Conifer defence against insects: microarray gene expression profiling of Sitka spruce (Picea sitchensis) induced by mechanical wounding or feeding by spruce budworms (Choristoneura occidentalis) or white pine weevils (Pissodes strobi) reveals large

Conifers are resistant to attack from a large number of potential herbivores or pathogens. Previous molecular and biochemical characterization of selected conifer defence systems support a model of multigenic, constitutive and induced defences that act on invading insects via physical, chemical, biochemical or ecological (multitrophic) mechanisms. However, the genomic foundation of the complex defence and resistance mechanisms of conifers is largely unknown. As part of a genomics strategy to characterize inducible defences and possible resistance mechanisms of conifers against insect herbivory, we developed a cDNA microarray building upon a new spruce (Picea spp.) expressed sequence tag resource. This first-generation spruce cDNA microarray contains 9720 cDNA elements representing c. 5500 unique genes. We used this array to monitor gene expression in Sitka spruce (Picea sitchensis) bark in response to herbivory by white pine weevils (Pissodes strobi, Curculionidae) or wounding, and in young shoot tips in response to western spruce budworm (Choristoneura occidentalis, Lepidopterae) feeding. Weevils are stem-boring insects that feed on phloem, while budworms are foliage feeding larvae that consume needles and young shoot tips. Both insect species and wounding treatment caused substantial changes of the host plant transcriptome detected in each case by differential gene expression of several thousand array elements at 1 or 2 d after the onset of treatment. Overall, there was considerable overlap among differentially expressed gene sets from these three stress treatments. Functional classification of the induced transcripts revealed genes with roles in general plant defence, octadecanoid and ethylene signalling, transport, secondary metabolism, and transcriptional regulation. Several genes involved in primary metabolic processes such as photosynthesis were down-regulated upon insect feeding or wounding, fitting with the concept of dynamic resource allocation in plant defence. Refined expression analysis using gene-specific primers and real-time PCR for selected transcripts was in agreement with microarray results for most genes tested. This study provides the first large-scale survey of insect-induced defence transcripts in a gymnosperm and provides a platform for functional investigation of plant-insect interactions in spruce. Induction of spruce genes of octadecanoid and ethylene signalling, terpenoid biosynthesis, and phenolic secondary metabolism are discussed in more detail.

Microarray expression profiling in melanoma reveals a BRAF mutation signature

We have used microarray gene expression profiling and machine learning to predict the presence of BRAF mutations in a panel of 61 melanoma cell lines. The BRAF gene was found to be mutated in 42 samples (69%) and intragenic mutations of the NRAS gene were detected in seven samples (11%). No cell line carried mutations of both genes. Using support vector machines, we have built a classifier that differentiates between melanoma cell lines based on BRAF mutation status. As few as 83 genes are able to discriminate between BRAF mutant and BRAF wild-type samples with clear separation observed using hierarchical clustering. Multidimensional scaling was used to visualize the relationship between a BRAF mutation signature and that of a generalized mitogen-activated protein kinase (MAPK) activation (either BRAF or NRAS mutation) in the context of the discriminating gene list. We observed that samples carrying NRAS mutations lie somewhere between those with or without BRAF mutations. These observations suggest that there are gene-specific mutation signals in addition to a common MAPK activation that result from the pleiotropic effects of either BRAF or NRAS on other signaling pathways, leading to measurably different transcriptional changes.

Saliva-derived DNA performs well in large-scale, high-density single-nucleotide polymorphism microarray studies

As of June 2009, 361 genome-wide association studies (GWAS) had been referenced by the HuGE database. GWAS require DNA from many thousands of individuals, relying on suitable DNA collections. We recently performed a multiple sclerosis (MS) GWAS where a substantial component of the cases (24%) had DNA derived from saliva. Genotyping was done on the Illumina genotyping platform using the Infinium Hap370CNV DUO microarray. Additionally, we genotyped 10 individuals in duplicate using both saliva- and blood-derived DNA. The performance of blood- versus saliva-derived DNA was compared using genotyping call rate, which reflects both the quantity and quality of genotyping per sample and the “GCScore,” an Illumina genotyping quality score, which is a measure of DNA quality. We also compared genotype calls and GCScores for the 10 sample pairs. Call rates were assessed for each sample individually. For the GWAS samples, we compared data according to source of DNA and center of origin. We observed high concordance in genotyping quality and quantity between the paired samples and minimal loss of quality and quantity of DNA in the saliva samples in the large GWAS sample, with the blood samples showing greater variation between centers of origin. This large data set highlights the usefulness of saliva DNA for genotyping, especially in high-density single-nucleotide polymorphism microarray studies such as GWAS.

Determining epithelial contribution to in vivo mesenchymal tumour expression signature using species-specific microarray profiling analysis of xenografts

Gene expression profiling using microarrays and xenograft transplants of human cancer cell lines are both popular tools to investigate human cancer. However, the undefined degree of cross hybridization between the mouse and human genomes hinders the use of microarrays to characterize gene expression of both the host and the cancer cell within the xenograft. Since an increasingly recognized aspect of cancer is the host response (or cancer-stroma interaction), we describe here a bioinformatic manipulation of the Affymetrix profiling that allows interrogation of the gene expression of both the mouse host and the human tumour. Evidence of microenvironmental regulation of epithelial mesenchymal transition of the tumour component in vivo is resolved against a background of mesenchymal gene expression. This tool could allow deeper insight to the mechanism of action of anti-cancer drugs, as typically novel drug efficacy is being tested in xenograft systems.

Meta-analysis of microarray data identifies GAS6 expression as an independent predictor of poor survival in ovarian cancer

Seeking new biomarkers for epithelial ovarian cancer, the fifth most common cause of death from all cancers in women and the leading cause of death from gynaecological malignancies, we performed a meta-analysis of three independent studies and compared the results in regard to clinicopathological parameters. This analysis revealed that GAS6 was highly expressed in ovarian cancer and therefore was selected as our candidate of choice. GAS6 encodes a secreted protein involved in physiological processes including cell proliferation, chemotaxis, and cell survival. We performed immunohistochemistry on various ovarian cancer tissues and found that GAS6 expression was elevated in tumour tissue samples compared to healthy control samples (P < 0.0001). In addition, GAS6 expression was also higher in tumours from patients with residual disease compared to those without. Our data propose GAS6 as an independent predictor of poor survival, suggesting GAS6, both on the mRNA and on the protein level, as a potential biomarker for ovarian cancer. In clinical practice, the staining of a tumour biopsy for GAS6 may be useful to assess cancer prognosis and/or to monitor disease progression.

Optimisation algorithms for microarray biclustering

In providing simultaneous information on expression profiles for thousands of genes, microarray technologies have, in recent years, been largely used to investigate mechanisms of gene expression. Clustering and classification of such data can, indeed, highlight patterns and provide insight on biological processes. A common approach is to consider genes and samples of microarray datasets as nodes in a bipartite graphs, where edges are weighted e.g. based on the expression levels. In this paper, using a previously-evaluated weighting scheme, we focus on search algorithms and evaluate, in the context of biclustering, several variations of Genetic Algorithms. We also introduce a new heuristic “Propagate”, which consists in recursively evaluating neighbour solutions with one more or one less active conditions. The results obtained on three well-known datasets show that, for a given weighting scheme,optimal or near-optimal solutions can be identified.

High-resolution SNP microarray investigation of copy number variations on chromosome 18 in a control cohort

Copy number variations (CNVs) as described in the healthy population are purported to contribute significantly to genetic heterogeneity. Recent studies have described CNVs using lymphoblastoid cell lines or by application of specifically developed algorithms to interrogate previously described data. However, the full extent of CNVs remains unclear. Using high-density SNP array, we have undertaken a comprehensive investigation of chromosome 18 for CNV discovery and characterisation of distribution and association with chromosome architecture. We identified 399 CNVs, of which loss represents 98%, 58% are less than 2.5 kb in size and 71% are intergenic. Intronic deletions account for the majority of copy number changes with gene involvement. Furthermore, one-third of CNVs do not have putative breakpoints within repetitive sequences. We conclude that replicative processes, mediated either by repetitive elements or microhomology, account for the majority of CNVs in the healthy population. Genomic instability involving the formation of a non-B structure is demonstrated in one region.

Rare chromosome disorders and their developmental consequences

Professionals working in disability services often encounter clients who have chromosome disorders such as Williams, Angelman or Down syndromes. As chromosome testing becomes increasingly sophisticated, however, more people are being diagnosed with very rare chromosome disorders that are identified not by a syndrome name, but rather by a description of the number, size and shape of their chromosomes (called the karyotype) or by a report of chromosome losses and gains detected through an advanced process known as microarray-based comparative genomic hybridisation (array CGH). For practitioners who work with individuals with rare chromosome disorders and their families, a basic level of knowledge about the evolving field of genetics, as well as specific knowledge about chromosome abnormalities, is essential since they must be able to demonstrate their knowledge and skills to clients (Simic & Turk, 2004). In addition, knowledge about the developmental consequences of various rare chromosome disorders is important for guiding prognoses, expectations, decisions and interventions. The current article provides information that aims to help practitioners work more effectively with this population. It begins by presenting essential information about chromosomes and their numerical and structural abnormalities and then considers the developmental consequences of rare chromosome disorders through a critical review of relevant literature.

Investigation of molecular mechanisms regulating biomineralization of pearl oyster Pinctada maxima

Biomineralization is a process encompassing all mineral containing tissues produced within an organism. The most dynamic example of this process is the formation of the mollusk shell, comprising a variety of crystal phases and microstructures. The organic component incorporated within the shell is said to dictate this remarkable architecture. Subsequently, for the past decade considerable research have been undertaken to identify and characterize the protein components involved in biomineralization. Despite these efforts the general understanding of the process remains ambiguous. This study employs a novel molecular approach to further the elucidation of the shell biomineralization. A microarray platform has been custom generated (PmaxArray 1.0) from the pearl oyster Pinctada maxima. PmaxArray 1.0 consists of 4992 expressed sequence tags (ESTs) originating from the mantle, an organ involved in shell formation. This microarray has been used as the primary tool for three separate investigations in an effort to associate transcriptional gene expression from P. maxima to the process of shell biomineralization. The first investigation analyzes the spatial expression of ESTs throughout the mantle organ. The mantle was dissected into five discrete regions and each analyzed for gene expression with PmaxArray 1.0. Over 2000 ESTs were differentially expressed among the tissue sections, identifying five major expression regions. Three of these regions have been proposed to have shell formation functions belonging to nacre, prismatic calcite and periostracum. The spatial gene expression map was confirmed by in situ hybridization, localizing a subset of ESTs from each expression region to the same mantle area. Comparative sequence analysis of ESTs expressed in the proposed shell formation regions with the BLAST tool, revealed a number of the transcripts were novel while others showed significant sequence similarities to previously characterized shell formation genes. The second investigation correlates temporal EST expression during P. maxima larval ontogeny with transitions in shell mineralization during the same period. A timeline documenting the morphologicat microstructural and mineralogical shell characteristics of P. maxima throughout larval ontogeny has been established. Three different shell types were noted based on the physical characters and termed, prodissoconch I, prodissoconch 11 and dissoconch. PmaxArray 1.0 analyzed ESTs expression of animals throughout the larval development of P. maxima, noting up-regulation of 359 ESTs in association with the shell transitions from prodissoconch 1 to prodissoconch 11 to dissoconch. Comparative sequence analysis of these ESTs indicates a number of the transcripts are novel as well as showing significant sequence similarities between ESTs and known shell matrix associated genes and proteins. These ESTs are discussed in relation to the shell characters associated with their temporal expression. The third investigation uses PmaxArray 1.0 to analyze gene expression in the mantle tissue of P. maxima specimens exposed to sub-lethal concentrations of a shell-deforming toxin, tributyltin (TBT). The shell specific effects of TBT are used in this investigation to interpret differential expression of ESTs with respect to shell formation functions. A lethal and sublethal TBT concentration range was established for P. maxima, noting a concentration of 50 ng L- 1 TBT as sub-lethal over a 21 day period. Mantle tissue from P. maxima animals treated with 50 ng L- 1 TBT was assessed for differential EST expression with untreated control animals. A total of 102 ESTs were identified as differentially expressed in association with TBT exposure, comparative sequence identities included an up-regulation of immunity and detoxification related genes and down-regulation of several shell matrix genes. A number of transcripts encoding novel peptides were additionally identified. The potential actions of these genes are discussed with reference to TBT toxicity and shell biomineralization. This thesis has used a microarray platform to analyze gene expression in spatial, temporal and toxicity investigations, revealing the involvement of numerous gene transcripts in specific shell formation functions. Investigation of thousands of transcripts simultaneously has provided a holistic interpretation of the organic components regulating shell biomineralization.